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.^ . Scaling level-spacing distribution functions in the "bulk of the spectrum" in ran- 

^^ . dom matrix models of N xN hermitian matrices and then going to the limit N ^ oo, 
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leads to the Fredholm determinant of the sine kernel sin7r(j; — y)/'ir{x — y). Simi- 
larly a double scaling limit at the "edge of the spectrum" leads to the Airy kernel 
_ . [Ai{x)Ai'{y) — Ai'(x)Ai(y)] /(x—y). We announce analogies for this Airy kernel of the 

,^ , following properties of the sine kernel: the completely integrable system of P.D.E.'s 

- ^-. , found by Jimbo, Miwa, Mori and Sato; the expression, in the case of a single interval, 

X, 

H , of the Fredholm determinant in terms of a Painleve transcendent; the existence of a 

commuting differential operator; and the fact that this operator can be used in the 
derivation of asymptotics, for general n, of the probability that an interval contains 
precisely n eigenvalues. 
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I. INTRODUCTION 

In this note we announce new results for the level spacing distribution functions obtained 
from scahng random matrix models oi N x N hermitian matrices at the edge of the sup- 
port of the (tree-level) eigenvalue densities when the parameters of the potential V are not 
"finely tuned." This universality class is already present in the Gaussian Unitary Ensemble. 
It is known [|T],0 that these distribution functions are expressible in terms of a Fredholm 
determinant of an integral operator K whose kernel involves Airy functions. 

There are striking analogies between the properties of this Airy kernel 

A{x)A'{y)-A'{x)Aiy) 

K[x,y) = 

x-y 

where A(x) = vA Ai(a;), and the sine kernel 

A sin7r(a; — y) 
TT X — y 

whose associated Fredholm determinant describes the classical level spacing distribution 
functions first studied by Wigner, Dyson, Mehta, and others |jl|. (In retrospect we should 
not have been surprised by this; the two kernels are, after all, both scaled limits of the same 
family of kernels.) We describe below three of these which we have found. The first is the 
analogue of the completely integrable system of P.D.E.'s of Jimbo, Miwa, Mori, and Sato 
when the underlying domain is a union of intervals. The second is the fact that in the case 
of the semi-infinite interval (s, oo) (the analogue of a single finite interval for the sine kernel) 
the Fredholm determinant is closely related to a Painleve transcendent of the second kind 
(the fifth transcendent arises for the sine kernel 0). And the third is the existence of a 
second order differential operator commuting with the Airy operator K. (The existence of 
such a differential operator in the sine kernel case has been known for some time 0,^.) This 
last fact leads to an explicit asymptotic formula, as the interval (s, oo) expands, for the 
probability that it contains precisely n eigenvalues (n = 1, 2, . . .) (the analogue of results in 



II. THE SYSTEM OF P.D.E.'S 
We set 

m 

-^ = U («2i-i, a2j) 
i=i 

and write -D(/; A) for the Fredholm determinant of K acting on I. We think of this as a 
function of a = (ai, . . . , a^m)- Then 

2m 

da\ogD{I-\) = -^(-iyi?(a„a,)c?a,- (1) 

i=i 

where R{x,y) is the kernel of the operator K{1 — K)~^. We introduce quantities 

q, = (1 - K)-'Aia,) , p, = (1 - K)-'A'ia,) , 

(which are the analogue of the quantities r^j of 0; see also 0) as well as two further 
quantities 

u= [A,{1-K)-^A) , v= [A,{1-K)-^A') 

where the inner products refer to the domain /. Then the equations read 

(7(3^ dj — Ct^j 

dak aj - ak 

aaj ^^. aj - ak 

— ^ = - ^^l-i)" ^ -^Pk + flj-gj + Pj« - ^qjv , 

COj ^_^^. aj - ak 

Moreover the quantities R{aj,aj) appearing in (|l|) are given by 

R{a,, a,) = Y.i-^r ^"^'^^ Z^''^'^' + P" - ^^^1 - 2P^-^^-" + ^^> ' 
k^j o-j a-k 

These equations are derived very much in the spirit of [Q; see also [^. 



III. THE O.D.E.'S 

For the special case I = (s, oo) the above equations can be used to show that g(s; A) (the 
quantity q of the last section corresponding to the end-point s) satisfies 

q" = sq + 2q\ {'= j-) (2) 

with g(s; A) ~ -\/A Ai(s) as s ^ oo. This equation is a special case of the P// differential 
equation P^TS], One can similarly derive for R{s) := R{s, s), which in view of (|I]) equals 

^ log D{I;X), 
as 

the third-order equation 

It is also the case that -R'(s) = —q{s; A)^ and this gives the following simple formula for 
D{r, A) in terms of a P// transcendent: 



D{r, A) = exp ( — / (x — s)q{x; A)^ j dx . 



This is much simpler than the corresponding representation of -D(/; A) for the sine kernel 
in terms of a Py transcendent. The fact that q{s; A) satisfies (|^) can also be obtained by 



combining some results in [^,0. Thus in this case of a semi- infinite interval our results 



have inverse scattering interpretations. 

IV. ASYMPTOTICS AND COMMUTING DIFFERENTIAL OPERATORS 

Again we take / = (s, oo) and consider asymptotics as s —>■ — oo. (Asymptotics as s — > oo 
can be obtained trivially from the Neumann series for (1 — K)~^.) From the random matrix 
point of view the interesting quantities are 






A=l 



This is the probabihty that exactly n eigenvalues lie in /. Using both differential equations 
(0) and (0), plus the fact that R' = — g^, we can obtain the asymptotics of i? as s ^ — oo: 

„, , 1 2 1 9 189 

^^^^^r -87 + 647-1287 + ---- 



(We also use the fact [0 that q{s; 1) is asymptotic to J—s/2 as s -^ — oo.) Therefore as 



-CX) 



where tq is an undetermined constant. (The analogue of this formula for the sine kernel was 
obtained by Dyson [|T^.) 

For asymptotics of E{n] s) for general n we introduce 

E{n;s) 

^^"'^^^=E(M- 

Successive differentiation of (|^) with respect to A, plus the known asymptotics of g(s; 1), 
allows us to find asymptotic expansions for the quantities 



qn[s) 



A=l 



9A" 

(for the analogue in the sine kernel case see 0); and these in turn can be used to find 
expansions for r{n;s). One drawback of this approach is that yet another undetermined 
constant factor enters the picture (in Toeplitz and Wiener-Hopf techniques, not available 
for the Airy kernel, fixed this constant). Another drawback is that when one expresses the 
r{n; s) in terms of the g„(s) a large amount of cancellation takes place, with the result that 
even the first-order asymptotics of r{n; s) are out of reach by this method when n is large. 
There is, however, another approach (briefly indicated in and with details in p] for 
the sine kernel case). We have 

-"-'v.J^.„ (i-A!:i:-:a-A.j « 

where Aq > Ai > ■ ■ ■ are the eigenvalues of the integral operator K (with A = 1). Now quite 
analogous to the fact that the operator with the sine kernel commutes with the differential 



operator for the prolate spheroidal wave functions ^0, is that the Airy operator commutes 
with the differential operator L given by 

Lf{x) = {{x - s) f\x))' - x{x - s)f{x) . 

An application of the WKB method, plus a trick, allows us to derive the following asymptotic 
formula for A,- with i fixed: 



1 _ A, ~ X^ 25^+3 ^=^V2+3/4 / _"^3/2 \ . ^ _2t ^ -OO). 



- - exp 

zl V 3 

(The analogue of this for the sine kernel is in [^.) From this it is seen that the term in @ 
corresponding to ii = 0, ii = 1, . . . , in = n — 1 dominates each of the others. In fact this 
term dominates the sum of all the others, and so 

Thus one can use (|^) to fix the constant in qi{s) mentioned above. 
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